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ABSTRACT 

In 2010 May 23-24, SDO observed the launch of two successive coronal mass ejections (CMEs), which 
were subsequently tracked by the SECCHI suite onboard STEREO. Using the COR2 coronagraphs 
and the heliospheric imagers (His), the initial direction of both CMEs is determined to be slightly 
west of the Sun-Earth line. We derive the CME kinematics, including the evolution of the CME 
expansion until 0.4 AU. We find that, during the interaction, the second CME decelerates from a 
speed above 500 km s _1 to 380 km s -1 , the speed of the leading edge of the first CME. STEREO 
observes a complex structure composed of two different bright tracks in HI2-A but only one bright 
track in HI2-B. In situ measurements from Wind show an "isolated" ICME, with the geometry of 
a flux rope preceded by a shock. Measurements in the sheath are consistent with draping around 
the transient. By combining remote-sensing and in situ measurements, we determine that this event 
shows a clear instance of deflection of two CMEs after their collision, and we estimate the deflection 
of the first CME to be about 10° towards the Sun-Earth line. The arrival time, arrival speed and 
radius at Earth of the first CME are best predicted from remote-sensing observations taken before 
the collision of the CMEs. Due to the over-expansion of the CME after the collision, there are few, 
if any, signs of interaction in in situ measurements. This study illustrates that complex interactions 
during the Sun-to-Earth propagation may not be revealed by in situ measurements alone. 
Subject headings: scattering — Sun: corona — Sun: coronal mass ejections (CMEs) 



1. INTRODUCTION 

The interaction of successive interplanetary coronal 
mass ejections (ICMEs) between the Sun and the Earth 
was first inferred from the analyses of multi-spacecraft 
in situ measurements during the 1980s, when the He- 
lios satellites in the inner heli osphere complemen ted the 
ISEE spacecraft near 1 AU flBurlaga et aL]|1987| >. With 
the improvement of coronagraphic observations, it was 
confirmed in the 1990s that successive ICMEs can merge, 
as well as interact with solar wind streams to form a 
compound stream or complex eject a ([G opalswamy et al. 



2001 1 |Burlaga et al.|2002| |20Q3| [Wang et al ~\M)2\ . Com- 



plex ejecta are often associated at E arth with extended 
periods of strong south ward B z (e.g. Wang et al.||2003[ 



Farrugia et al.| |2006a| and with intense geomagnetic 
storms (Burlaga e t al.| 19871 Farrugia et al.l2006a|b| |Xie| 
et al.|200 6|). Until the launch ol the Solar- 1 errestrial Re- 



lations Observatory (STEREO, see: |Kaiser et al.|2008[ ) m 



2006, there was very limited data available beyond the 
field-of-view (FOV) of coronagraphs (~ 0.15 AU). In ad- 
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dition, in situ measurements outside of Earth's direct 
vicinity have been limited to planetary missions, which, 
generally, cannot easily be combined with measurements 
at 1 AU to study the interaction of successive CMEs 
and the formation of complex ejecta. Most of what we 
have learned about CME-CME interaction comes from 
nu merical simulat i ons. S uch simulations were pioneered 
by jVandas et aT] fll997| and [Schmidt fc Cargill| ( |2004 ) 
and they have been performed si nce using 2. 5-D and 3-D 



magneto- hydro dynamica l code s (|Wu et al.| 2002 ; Odstrcil 
|et al.||2003| |Lugaz et al]|2005| pGong et al.||2U(j|) |2U 



A lew real events have also been simulated ( |Wu et a. 
2007||Lugaz et al.|2007] ). 

Based on these simulations, it was found that the 
interaction of successive CMEs can cause intervals of 
southward B z with longer duration and larger magnitude 



flWang et al.|2003[ Lugaz et al.|2008a |Xiong et al.|2 009 ) 

and a lso increase the density in the sneath ( |Lugaz et al.| 



2005 ), which could have important c onsequences lor geo- 
etlectiveness (Farrugia et al. 2006a). Some of the main 



remaining areas ol investigation related to CME-CME in- 
teraction are: (i) the momentum exchange between the 
successive eruptions during their interaction, (ii) the fate 
of the related shocks, (hi) the possibility of CME-CME 
"cannibalism" through total reconnection of one of the 
CMEs, and, (iv) the deflection and rotation of a CME 
during the interaction. Regarding the momentum ex- 
change during the interaction, it has been proposed that 
the interaction is elastic, perfectly inelastic, slowe r than 
inelastic due to the presence of the MHD obstac le ( |Tem-| 
mer et al.|2012[ ), super-elastic (Sh en et al.|2012[ ), or that 



the momentum is transferred through the propaga tion of 
the shock insi d e the fir st CME (see dis cussi ons in|Lugaz 



et ah] (|2009aJ), |Farrugia et al.| ( |2006a[ ) and |Xiong et al. 
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Fig. 1. CMEs of 2010 May 23 (CME1, top) and 2010 May 24 (CME2, bottom) as observed by SDO 193 A (left), COR-2 onboard 
STEREO- A (middle) and STEREO-B (right). 



(2006)). 



With the Sun-Earth Connection Coronal and He- 



liospe ric Investigation (SECCHI) suite (Howard et al 
2008) onboard STEREO, it is now possible to track the 
density fronts a ssociated with CM Es from the corona to 
1 AU (e.g., see Davies et al.|2009|, to fo recast the ar- 
rival at Earth or the 1CME ( |Davis et al.||2011[ ). A par- 
ticularly promising line of research is to combine these 
remote-sensing observations with in situ measurements 
from ACE, Wind, and the S TEREO spacecraft to an- 
alyze the evolution of CMEs ([Wood et al 1120091 |Mostl | 



et al 120101 [Liu et al.|2010bl |Romllard|^(ji(j| |Rollett etaT 



2012). Because of the deep solar minimum encompass 
mg most of the STEREO mission, there have not been 
many instances of successive and potentially interacting 
CMEs observed by SECCHI until early 2010. One ma- 
jor exception is the 2007 January 24-25 CMEs. This 
event has been fully analyzed in a number of articles 
using STEREO and SMEI data ([Harrison et aT] [2009 



Webb et al.||2009l |Lugaz et al . 2009a) as we l l as numer- 
ical simulations ([Eugaz et al.||2008b[ |2009b[ |Odstrcil fc 



Pizzo 20091). While this was an instance of interacting 
CMEs, the lack of in situ data and a long data gap make 
this event not particularly enlightening to understand 
CME-CME interaction. Another example of recent ob- 
servations of multiple CMEs is the series of eruptions 
in late July and early August 2010. These CMEs have 
been extensively studied ([Harrison et al.|[2012) |Liu et al. 

2012p, but the number of involved 



2012 
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CMEs (up to six, tour ol which occurred in one day) and 
the complexity of the observations make these events not 



optimal for studying the physical processes underlying 
CME-CME interaction. 

In this paper, we present an analysis of observa- 
tions by STEREO/SECCHI combined with measure- 
ments from the Wind spacecraft of two CMEs, which 
erupted on 2010 March 23 and 24. In section [2j we dis- 
cuss the remote-sensing observations of the two CMEs 
in extreme-ultraviolet (EUV) i magers, coro nagraphs and 
heliospheric imagers (His, see: Eyl es et al.||2009| ). In sec- 
tion [3j we derive the CME directions using a number of 
methods. We analyze the measurements and derive the 
CME kinematics in the His FOVs in section [4j and we 
analyze the associated in situ measurements in section 
[5| We discuss our results and propose a coherent sce- 
nario for the CME-CME interaction in section [6] and we 
conclude in section 

2. REMOTE-SENSING OBSERVATIONS 

2.1. First CME 

On 2010 May 23, SDO/AIA and SOHO/EIT observed 
a Bl flare and a filament eruption originating from 
N19W12. STEREO-A was at a heliocentric distance of 
0.956 AU and an angular separation of 71.5° to the west 
of Earth. STEREO-B was at a distance of 1.014 AU and 
an angular separation of 70° east of Earth. An erup- 
tion (hereafter, CME1) was detected first by COR1-B at 
16:05 UT where it appeared as a western limb eruption 
(from the spacecraft perspective). It was detected by 
COR1-A at 16:35 UT as a eastern limb eruption. Ob- 
servations of CME1 in SDO, COR2-A and COR2-B are 
shown in the top row of Figure 1. The CME appears to 



A COMPLEX EVENT OBSERVED BY SECCHI 







Fig. 2. Background-subtracted images of the 2 CMEs on 2010 May 25 in HI1 onboard STEREO-A (left) and STEREO-B (right). The 
two rows correspond to 00:09 UT (top) and 06:49 UT (bottom). The concentric circles show the elongation angles and the radial lines the 
position angles, both in degrees. The position in elongation angle of CME1 is around 17° (top) and 21° (bottom); that of CME2 is around 
10° (top) and 15° (bottom). [This figure is available as two mpeg files in the electronic edition of the Journal]. 



be deflected slightly southward towards the ecliptic as it 
propagates through the COR1 FOV. It entered the FOV 
of COR2-A at 17:24 UT and of COR2-B at 17:54 UT. 
The radial speed in COR2-A was 362 km s _1 and the 
central position angle (PA) about 95°. Speed and cen- 
tral PA in COR-2 quoted here and thereafter are from 
the CACTUS database ( |Robbrecht fc Berghmans||2004| ) 
unless noted otherwise, and the speeds are plane-oi-the- 
sky speeds. In COR2-B, the speed was 378 km s _1 and 
the central PA about 275°. CME1 was also observed by 
LASCO/C2 starting at 18:30 UT as a faint asymmetric 
halo with a central PA of 320°. The average speed of 
the halo CME was 278 km s _1 (values from the CDAW 
CME list). We det ermined the mass of this CM E us- 
ing the procedure of Colaninno & Vourlidas (2009). We 



found a mass of 1.5±0.1 x 10 id _ 

CME1 first entered the HI1-A FOV at 21:29 UT on 
May 23 and HI1-B at 22:09 UT. It had is a typical 3- 
part structure with a dense core that entered the HI1-A 
and HI1-B FOVs at 02:09 UT on May 24. The CME first 
entered HI2-A a t 22:09 UT on May 24 and HI2-B at 00:09 



UT on May 25. |Song et al (2012) recently reported and 
analyzed a series oi blobs propagating and interacting in 



the corona within the current sheet behind CME1 from 
04:00 to 13:00 UT on May 24, as well as an associated 
type-Ill radio burst detected by STEREO-A at around 
10:24 UT on May 24. There was also an unrelated narrow 
CME at around 06:00UT on May 24 observed at PA 66° 
in COR2-A and at PA 310° in COR2-B. 

2.2. Second CME 

The second CME (hereafter, CME2) studied here was 
associated with a Bl flare and a filament eruption from 
N18W26 observed by SDO/AIA and SOHO/EIT. It was 
the same filament channel erupting as for the previous 
CME, but due to solar rotation it had moved by about 
12° more to the west. An eruption was detected first by 
COR1-B at 13:05 UT on May 24, where it appeared as 
a western limb eruption. It was detected by COR1-A 
at 13:45 UT as an eastern limb eruption. Observations 
of CME2 in SDO, COR2-A and COR2-B are shown in 
the bottom row of Figure 1. CME2 entered the FOV of 
COR2-A at 14:24 UT and the FOV of COR2-B at 14:54 
UT. The speed in COR2-A was about 500 km s -1 and 
the central position angle (PA) about 95°. In COR2-B, 
the speed was about 520 km s _1 , and the central PA ap- 
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Fig. 3. Graduated Cylindrical Shell (GCS) fitting of CME1 (top) and CME2 (bottom). The images are from STEREO-B /COR2, 
LASCO/C2 and STEREO-A/COR2, from left to right and the green overlay shows the fitted shape. See texts for details. 



proximately 270°. These speeds correspond to the "core" 
of the CME, which propagated to the south of the eclip- 
tic (near PA 110°, as seen from A and PA 250°, as seen 
from B). In addition to this structure, which was rela- 
tively narrow (40° wide in PA), there was a much wider 
(90° wide) and faster "leading edge" centered around the 
ecliptic whose speed was around 650 km s" 1 . CME2 was 
also observed by LASCO/C2 starting at 14:06 UT on 
May 24 as a faint asymmetric halo with a central PA 
around 280°. The average speed o f the halo CME was 
427 km s - 1 . Using the procedure of Colaninno fc Vourli- 



daslfl2009| ), we find a mass of 1.0±0.1 x 10 iD g for CME2. 
CME2 entered the HI1-A FOV at 18:09 UT on May 24 
and HI1-B at 17:29 UT. It entered the HI2-A FOV at 
06:09 UT on May 25 during a data gap of STEREO- 
B/SECCHL 

2.3. Indications of a Probable Interaction 

Assuming constant speeds of 380 km s _1 and 
500 km s- 1 (as reported on the CACTUS CME catalog) 
for the two CMEs, the two leading edges should cross 
at around 0.8 AU. Assuming speeds of 400 km s _1 and 
650 km s _1 (speed of the fastest fronts), the interaction 
should have happened by 0.55 AU. While a deceleration 
of the CMEs is expected, interaction between the front 
of CME2 and the back of CME1 is anticipated to hap- 



pen much closer to the Sun than the distance estimated 
above. We also note that the two CME fronts overlap 
in the HI2 FOV, possibly indicating interaction of the 
CMEs. Background-subtracted images of the two CMEs 
in the HI1 FOV are shown in Figure 2, corresponding 
to 2010 May 25 at 00:09 UT (top row) and 06:49 UT 
(bottom row). 

3. DIRECTION OF THE BRIGHT FRONTS 

Determining the directions of the bright fronts ob- 
served by SEC CHI is essential for understanding the 
kinematics and evolution of the two CMEs. Southward 
deflection away from the flare site was observed for both 
CMEs. Therefore, we cannot simply use the flare loca- 
tion as the direction of propagatio n of the CMEs. Single- 
spacecraft fitting metho ds (Sheeley et al. 1 999} [Rouillard| 



et al. 



2008 



Lugaz||2010| IDavies et aT| 2012) can only be 



used ii the speeds of the CMEs are relatively constant. In 
the case of a CME-CME interaction, there would be large 
variations in the CME speeds. Finally, it is not even clear 
which bright front in HI2 corresponds to which CME. 

For these reas ons, we rely on stereoscopic methods: 
triangulation (see Liu et al.|2010a[ ) , which is based on the 
Fix ed-<E> (F3>) geometry ; and tangent-to-a-sphere method 
(see |Lugaz et al.pOlO ), which is based on the Harmonic 
Mean (HM) geometry. Triangulation has been mostly 
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Fig. 4. — J-maps for the two CMEs and the discussed tracks. Left: STEREO-A, right: STEREO-B. Top: Running-difference, bottom: 
background-subtracted. From left to right, there is the leading edge of CME (shown in blue), its black edge (BE1, red), the end of the 
leading edge or beginning of the cavity (orange), the end of the CME1 (pink), the leading edge (green) and black edge (yellow) of CME2 
(in red) as shown in the J-maps. Brown is used to highlight a secondary front in HI2-A discussed in the text. 



used in the EUVI, COR1, COR2 and HI1 FOV s, and it 
has been shown (e.g., see jde Koning et al.||2QQ9[ ) to yield 
reliable results as compared to other methods, .because it 
assumes that the same plasma element is being observed 
by both spacecraft, triangulation is expected to break 
down at large elongation angles and for wide CMEs. The 
tangent-to-a-sphere method takes into account the fact 
that the same part of the CME is not imaged by the 
two spacecraft. It is expected to yield better results for 
those CMEs observed as a halo by one of the spacecraft, 
or for which the width cannot be neglected. We also 
use the visual fitting of the CMEs in the coronagraphs 
FOV using the Graduated Cylindrical Sh ell (GCS) model 
( jThernisien et al.|2QQ9| |Thernisien|2Q 1 1 [ ) , which assumes 
a Hux rope shape. In the next section, we analyze the 
observations in the COR2 and the HI1 FOVs to derive 
the CME directions. 

3.1. Direction of the CMEs Using the GCS Model 

We use the three coronagraphic views on CME1 at 
around 20:30 UT on May 23 from STEREO-A/COR2 
(at 20:24 UT), STEREO-B/COR2 (at 20:24 UT) and 
LASCO/C2 (at 20:30 UT). The three views are shown 
in the top row of Figure 3. At this time, the CME was 
at a distance of around 9.5 Rq. The best visual fit for 
CME1 is a direction of S0W10 with a large tilt angle of 



65° with respect to the ecliptic, and a small half-angle of 
15°. For a disc ussion of the unce r tainti es of the method, 
please refer to |Thernisien et al.] (|2009| ) . Typical uncer- 
tainties are in the range of ±5 U , =tlU u , and ±20° for the 
longitude, the half-angle, and the tilt, respectively. 

For CME2, we use the three views at around 16:30 UT 
on May 24, at which time CME2 was also around 9.5 R@ 
(the three views are shown in the bottom row of Figure 3, 
corresponding to 16:24 UT for the COR2 views and 16:30 
UT for the LASCO view). The best fit for CME2 is a 
direction of S02W26, with a tilt angle of 60° and a half- 
angle of 20°. As can be seen from Figure 3, for both 
CMEs, the LASCO view is essential for constraining the 
direction of the CMEs since both CMEs had a similar 
aspect in STEREO-A and STEREO-B. 

3.2. Direction of CME1 in Ell FOV 

In the HI1 FOV, we first track three features associated 
with CME1: the leading edge (as seen in background- 
subtracted or running-difference J-maps); the "black" 
edge (the transition from white to black, see below) 
as seen in running-difference J-maps; and the back of 
the leading edge as observed in background-subtracted 
J-maps. J-maps are time-elongati on maps created fol- 
lowi ng the procedure ex plained in |Sheeley et al. (1999) 
and Davies et al. (2009). In the bottom panels of Fig- 
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Fig. 5. — Direction of the first (top) and second (middle) CMEs. 
Crosses are used for results from the triangulation method and 
diamonds for the tangent-to-a-sphere method. The leading edges 
are shown in red, the "black" edges in black and the back of the 
CME in blue. The bottom panel shows a sketch of the geometry of 
the observations on 2010 May 24-27 as seen from the ecliptic north. 
The initial directions (12° and 28°) are shown with solid arrows 
and the deflection during the interaction is shown with dashed 
arrows (see text for details). 



ure 4 (background-subtracted J-maps), the leading edge 
of CME1 is highlighted in blue and the back of the 
leading edge in orange. In the top panels of Figure 4 
(running-difference J-maps), the "black" edge is high- 
lighted in red. 

The "black" edge corresponds to the apparent end of 
the leading edge as seen in the running-difference images. 
This is easier to track because the contrast is optimal 



between bright and dark regions as compared to the real 
leading edge, which is between grey and bright regions. 
However, it is not clear to what physical part of the CME 
front (if any) the "black" edge corresponds. By tracking 
three different features, we attempt to determine if all 
parts of a CME can be reasonably assumed to have the 
same direction. 

Using triangulation, the propagation direction of the 
leading edge is derived to be about 14° (positive angles 
refer to west of the Sun-Earth line), the direction of the 
"black" edge is about 10.5° and that of the back of the 
CME is about 10.5°. In all cases, the direction remains 
almost constant with distance until about 14-15° elonga- 



tion w hen it increases. Using the method of |Lugaz et al. 
pOlo] ), the CME direction is found to be 27.5°, 20.5° 
and 20.5° for the leading edge, the "black" edge, and the 
back of the CME, respectively. Results using HI1 data 
are summarized in the top panel of Figure 5. 

This Figure shows two important results. 1) The direc- 
tion as derived using the actual leading edge, the "black" 
edge, or the back of the CME are within the typical 
error bar of the methods (±5°), with the direction ob- 
tained from the "black" edge being more steady than 
that derived from the actual leading edge. This validates 
all previous studies which have typically used the black 
edge of the CME to derive its direction and kinematics. 
2) Our study shows that, for the first 24 hours of its 
propagation in the heliosphere (past 20 i?©), the CME 
has a constant direction, independently of the method 
used to derive this direction. Finally, it should be noted 
that 3) the CME heliocentric distance in the HI1 FOV 
has only a very weak dependence on the method cho- 
sen to derive the propagation direction. This has been 
pointed out before, especially for fitting methods; both 
the HM and F<1> fitting tend to produce different direc- 
tions and speeds, but s imilar arrival times (for e x ample 



see 



H oward fc Tappin (2009), Mostl et al. (2011), Lugaz 
eTal.| ( |201^| )).~ 



We believe that the direction of CME1 derived by tri- 
angulation should be trusted because: i) it corresponds 
to the approximate location of the source region of the 
CME (W12) as well as the direction obtained from the 
GCS fitting, and, ii) the CME i s relatively narrow and 
the assumptions of the method of Lugaz et al. (2010) are 
not valid in this case. For the rest oi the analysis, we 
use a direction of 12° (the average of the direction from 
triangulation of the three tracked features) and we use 
the F<£ approximation in HI1 to determine the distance 
and kinematics for all tracks associated with CME1. 

3.3. Direction of CME2 in Ell FOV 

We follow the same procedure as described above for 
CME2. It should be noted that the edges of CME2 are 
less well defined because of the presence of bright features 
associated with CME1. In the bottom panels of Figure 4, 
the leading edge of CME2 is highlighted in green and, in 
the top panels of Figure 4, the "black" edge of CME2 is 
highlighted in yellow. Using triangulation, the direction 
of the leading edge is about 15° and the direction of the 
"black" edge is about 14°. Using the method of [Lugaz] 
et al.| fl2010| ), the CME direction is 29° and 27° for the 
leading and the "black" edges, respectively. 

We believe that the direction of the CME derived by 
the tangent-to-a-sphere method should be used, because 
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it corresponds to the approximate location of the source 
region of the CME (W26) as well as the direction ob- 
tained from the GCS fitting (W26). Hereafter, we use the 
HM approximation to analyze observations of CME2 in 
HI1. A sketch of the geometry including the two CMEs, 
as seen from ecliptic north, is shown in the bottom panel 
of Figure 5. 

4. EVOLUTION OF THE CMES IN HI1 AND HI2 FOV 
4.1. HI Data 

In addition to the leading edges (real and black, as 
discussed above), we track the back of the leading edge 
of CME1 as seen in background-subtracted images (or- 
ange track in Figure 4). At 1 AU, magnetic ejecta typ- 
ically have low density ( jBurlaga et al. 1981). Magnetic 
flux ropes are thought to be present in t he dark cavi 



ties o f three-part CME structures (e.g., see Hundhausen 
|1993[ ). Ejected high-density filament material (the CME 
core) is at first embedded inside this cavity and, at later 
times, it is expected to over-expand and become part of 
the low-density region. Therefore, the back of the lead- 
ing edge can be assumed to correspond appro ximately 
to the front of the magnetic ejecta. Recently, Howard | 



fc DeForest (2012) have shown how a magnetic cloud 
observed at Earth can be traced back to the dark cav- 
ity as seen in background-subtracted images in corona- 
graphs and heliospheric imagers. Density enhancements 
or V-shape structures are often observed by SMEI and 
the His corresponding to the back of the magnetic ejecta 
dKahler fc Webb||20Q7| |Harrison et"al|l2QQ8[ ). We also 
track one such feature at the back of ClvlEl (pink track 
in Figure 4), which corresponds to the compression at the 



back of t he cavity or magnetic ejecta ( |Howard fc DeFor- 
est| 2012 ). The measurements of the rive tracked fronts 



m Hll are summarized in the top left panel of Figure 6. 

We find that the temporal offset between the leading 
edge (asterisks) and the black edge (triangles) increases 
from 2 hours to about 3.3 hours as CME1 propagates in 
the HI1-A FOV. In the HI2-A FOV, we track the ap- 
parent brightness maximum (or center of the track) in 
the running-difference J-map as well as the black edge, 
and we find a temporal offset between the two tracks of 
about 3 hours. Assuming that the apparent center of 
the track is the actual center of the leading edge, this 
means that the temporal offset between the leading edge 
and the black edge of CME1 is about 6 hours. An off- 



set of about 4-6 hours was also estimated in Lugaz et al. 
(2012) and, in general, it should be corrected in studies 
of CME arrival times using the black track. The offset is 
approximately equal in HI1-B to what is found in HI1-A. 

This finding means that the black edge seen in running- 
difference images and J-maps is not simply a result of the 
running-difference procedure, since the offset between 
the CME leading edge and the black edge is not constant 
within an instrument. The increasing offset is most likely 
due to the expansion of the CME sheath. However, there 
is a strong difference in the value of the offset between 
Hll (~ 2 — 3 hours) and HI2 (~ 4 — 6 hours), which in- 
dicates that this offset is also partially controlled by the 
resolution and frequency of the observations. 

4.2. Determining the distances of the CMEs in the 
inner heliosphere 



In the top right panel of Figure 6, we show the time- 
distance evo l ution of the different tracks. As done in 
|Lugaz et al. (2010), we determine the typical error as- 
sociated witn the manual selection of the front by esti- 
mating the accuracy of the elongation measurements (as 
shown in the top left panel of Figure 6) to be ±0.15° 
(only valid for Hll data). Then, the error on the dis- 
tances is typically ±0.6 Rq and the error on the velocity 
is ±85 km s _1 . In the bottom left panel of Figure 6, 
we show the speed of the leading edge of the two CMEs; 
we use the black edge of CME1, because the results are 
steadier, and the leading edge of CME2 because the black 
edge is not well defined early on (see the yellow track in 
Figure 4). 



4.3. Width of the First CME 

Using the data obtained from HI1-A, we calculate some 
of the radial widths (or radial extents) associated with 
CME1: total ICME width, total ejecta width (without 
the sheath), and sheath width. Based on the typical 
error for the distance, the widths have an uncertainty 
of ±1.2 R(?>. A sim ilar study was recently completed by 



Savani et al.| ( |2Q12| ) for four CMEs observed by SECCHI. 
As can be seen m the bottom right panel of Figure 6, 
the three widths increase with distance, except that the 
ejecta and total ICME width become constant or start to 
decrease on May 25 at 00:09UT. This does not appear to 
be connected with the front of the CME, as is clear from 
the fact that the width of the sheath region continues 
to increase. This marks the probable beginning of the 
interaction between the two CMEs. 

We fit the widths, W as power laws with respect 
to the radial distance, r before the start of the CME- 
CME interaction. We find the following relationships: 

W lC ME = 0.30r° c 7 6 E , W ejecta = 0.21r° e 8 c 2 ta and Sheath = 

0.11rjQ 5 j^ E , where W and r are both expressed in AU. The 
value of the power index for the width of th e ejecta and 
the ICME are very close to the relation of Both mer fc| 
|Schwenn| ( |1994[ ), who found a power law of 0.78 based on 
in situ measurements at 1 AU for magnetic clouds, and 
relatively close to the powe r-law values of 0.92 and 0.9 1 



found by |Liu et al.| ( |2005D and |Gulisano etaL| |2010j) 
respectively! All these studies also find a typical value 
for the ejecta width at 1 AU of 0.24±0.02 AU, close 
to the value of 0.21 AU found here. Here, contrary to 
many other studies, the sheath is not made of shock- 
compressed solar wind, because CME1 is probably too 
slow to drive a shock. What is referred to as the sheath 
is a combination of swept-up solar wind mass and possi- 
bly some coronal mass (part of what initially comprises 
CME1). 

Finally, we make a crude estimate of the angular half- 
width of CME1 in ecliptic longitude. To do this, we 
assume that the radial width as derived here is equal 
to the CME longitudinal width. This assumption is not 
arbitrary, since, in fact, the axis of CME1 is found to 
be highly tilted using the GCS fitting (see also measure- 
ments at 1 AU) and CME cross-sections close to the Sun 
are expected to be nearly circular. Under these assump- 
tions, we find that the angular half-width of CME1 is 
about 8-10°. A different, and more rigorous, estimate 
ca n be obtained using the stereoscopic method (model 2) 
of Lugaz et al. ( 2010| ), which allows to derive the actual 
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Fig. 6. — Top left: Time-elongation measurements in HI1-A of the tracked fronts, corresponding to the leading edge of CME1, its black 
edge, the end of the leading edge or beginning of the cavity, the end of the CME1, all in black, and the leading edge and black edge of 
CME2 (in red) as shown in the J-maps. Top right: Heliocentric distances of the tracked fronts for CME1 (black) and CME2 (red). The 
fronts are the same as described in top left with the same conventions. Bottom left: Speed of the black edge of CME1 and leading edge of 
CME2. Bottom right: Widths of CME1; diamonds: total ICME (sheath ± ejecta), stars: sheath, and square: ejecta. Error estimates are 
±0.6 Rq for the distances, ±1.2 Rq for the width, and ±85 km s _1 for the velocities (see text for details). 



CME1 longitudinal half- width under a set of assump- 
tions. Using HI1 data and a fixed direction of 12°, we 
find that the half-width is 11° ±2°, in relatively good 
agreement with what we derived above. Such a small 
half-width is consistent with the use of the F<3> approx- 
imation for this CME and has some important conse- 
quences when we consider the physical processes during 
the CME-CME interaction (section 6). 

4.4. Kinematics of the CMEs in Ell FOV 

Initially on May 24, the leading edge of CME1 prop- 
agates with a speed of about 360 km s _1 and the lead- 
ing edge of CME2 with a speed of 620 km s -1 , corre- 
sponding approximately to the speed derived in COR2. 
A direct consequence of the difference in speeds is that 
CME2 quickly catches up with CME1. In addition to 
the speeds plotted in Figure 6, we find that the begin- 
ning of the magnetic ejecta (cavity) of CME1 has a speed 
of 340 km s _1 and that the end of CME1 has a speed of 
about 260 km s -1 , reflecting the expansion of the mag- 
netic ejecta (as described in the previous section). 

At 00:09 UT on May 25, CME1 starts to contract in 
radial extent, with the front and rear of CME1 becom- 
ing closer together, and at around this time (01:19 UT 



on May 25), the leading edge of CME2 decelerates to 
reach a near constant value of 380 km s _1 at 02:40 UT. 
This supports our claim that, indeed, the leading edge of 
CME2 has collided with the back of the magnetic ejecta 
of CME1. The material that makes up the leading edge 
of CME2 cannot propagate through the magnetic bar- 
rier and it is forced to propagate with the same speed 
as CME1. New material is not expected to accumulate 
between the two CMEs but the existing material cannot 
easily be removed until the two CMEs separate, if they 
do. The speed of this dense material is now constrained 
by the speed of CME1. 

4.5. Evolution of the CMEs in the HI2 FOV 

The change in temporal and spatial resolution as 
well as intensity threshold and the difficulty to use 
background- subtracted images, make it complicated to 
follow the evolution of the CMEs in the HI2 FOV. In 
addition, the long data gap in HI-B further muddles the 
picture. From a visual inspection of the J-maps, the first 
striking feature in HI2 is the fact that there are two bright 
tracks associated with the CMEs in HI2-A but only one 
bright track in HI2-B. Below, we quickly discuss the ori- 
gin of these tracks. 
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Fig. 7. — STEREO-A HI2 running-difference images on May 25 showing the "splitting" of the first front and its "interaction" with the 
second front and also illustrating the fact that it appears that after on May 26, CME2 is now ahead of CME1. 



The idea that the first and second tracks in HI2-A cor- 
respond to different CMEs before and after the short data 
gap in HI-A on May 26 can be inferred by inspecting the 
actual images (instead of the J-map). Throughout its 
propagation, CME1 appears relatively symmetric in po- 
sition angle with respect to the ecliptic, whereas CME2 
has an ear-shape with a flattening or change in indenta- 
tion around position angle (PA) 110° (see Figure 2, for 
example). As shown in Figure 7, on May 26, the first 
CME observed in HI2-A has the shape characteristics of 
CME2. 

We want to emphasize the following: while running- 
difference images in HI2-A show two (or three) bright 
tracks at all times, separated by dark regions, until the 
beginning of May 26, there is in fact only one wide bright 
front as seen in background-subtracted images (as can 
be seen on the bottom left panel of Figure 4). In HI2-A 
at 14:09 UT on May 25, there is a "split" of the first 
track visible at around 25° (see left and middle panels of 
Figure 7). It appears as if part of the first track becomes 
part of the second track (this transient third track is 
highlighted in brown in the running-difference J-maps). 
Finally, it is worth noting that the first track in HI2- 
A clearly has some geometr ical acceleration ( simila r to 
what has been discussed in iRouillard et al.| ( |2QQ8[ ) for 
example). This, combined with the fact that no similar 
front is observed in STEREO-B images, shows that this 
front propagates in a direction away from the Sun-Earth 
line, more towards STEREO-A, as CME2 is sketched in 
the bottom panel of Figure 5. 

Overall, images and J-maps from HI2 strongly sug- 
gests that the first track in HI2-A corresponds to CME2, 
which is at a larger heliocentric distance and closer to 
STEREO-A than CME1. We can hypothesize that the 
"split" of the first track in HI2-A, which merges with the 
second track, occurs at the time when the elongation an- 
gle of CME2 becomes equal and then larger than that of 
CME1. 

5. ARRIVAL TIME AT 1 AU AND IN SITU 
MEASUREMENTS 

5.1. Overview of the In Situ Measurements 

Magnetic field and plasma data from the Wind space- 
craft are shown in the left column of Figure 8 for the 
period 20 UT, May 27 to 8 UT, May 30, 2010. Wind was 



orbiting the LI Lagrangian point and at the start and end 
of the period it was at (262.2, -13.0, 25.3) R E and (261.9, 
-18.2, 25.2) R E , respectively (GSE coordinates). The 
plasma data are from th e Three-Dimensional Plasma An- 
alyzer (Lin et al. 1995), and the magnetic field dat a are 



from the Magnet ic Field Investigation instrument ( |Lep-| 
|ping et aL [1995). Both have a temporal resolution of 
three seconds! The panels show, from top to bottom, 
the proton plasma density (in red, the alpha particle- 
to-proton number density ratio in percent), the proton 
temperature (in red, the expected p roton temper ature 
for normal solar wind expansion after |Lopez] (|1987 ) ) , the 
bulk flow speed and the east- west component (y y ) of 
flow vector, the pressures (total, P t , magnetic, P m , and 
the sum of the proton and electron thermal pressure, 
P p + P e ), the strength and GSE components of the mag- 
netic field, B, colored as shown, and the proton plasma 

The data reveal a clear, single interplanetary coronal 
mass ejection passing Wind at the times shown between 
the second and third vertical guidelines. It is character- 
ized by a strong field, smooth rotations of B components, 
and a low proton /3. The alpha particle-to-proton density 
ratio is generally also higher than average. The proton 
temperature is comparable and even above the expected 
proton temperat ure, so the config uration is not a mag- 



netic cloud (Burlaga et al. 1981). The rotation of the 



field indicates, however, that it is a magnetic flux rope 
(see further below). 

A shock (marked as S) precedes the ICME by ~18 
hrs. We examined this shock by th e velocity coplanarity 
method ( Abraham- Shrauner| 19 72 ) and obtained a shock 
speed of 392 km s _i . This is comparable with the speed 
at the leading edge of the ICME, consistent with it being 
driven by the transient. The angle between the upstream 
magnetic field and the shock normal is 14.2°, hence it is 
a quasi-parallel shock. While it may be surprising that 
such a slow magnetic cloud drives a shock, it should be 
noted that the Alfvenic and sonic speeds upstream of 
the shock are both below 40 km s _1 so that the CME 
is faster than the fast magnetosonic speed in the solar 
wind frame. The sheath preceding the CME is found to 
be about 0.16 AU long. 

We examined the ICME magne tic field data by a 
minimum variance technique (e.g. Sonnerup & Cahill 
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Fig. 8. — In situ data as observed by the Wind spacecraft (left) from May 27 20 UT to May 30 08 UT and zoom-in on the ejecta and 
the last part of the sheath after minimum- variance analysis (top right). The bottom right panel shows the Grad-Shafranov reconstructed 
cross-section map (see text for details). The black contours show magnetic field lines in a plane perpendicular to the flux rope axis, which 
is indicated by a white dot; and the axial field strength is color coded. Along the Wind trajectory, remaining velocities in the deHoffmann- 
Teller frame are indicated by green arrows, and observed magnetic field components in the plane of the cross-section are shown with white 
arrows. 



1967) to confirm the visual inspection. With a ratio of Gr ad- Shafr anoy reconstruction ( Hu fc Sonnerup| 2001 



intermediate-to- minimum eigenvalues of 8.1, the routine 
returns a reliable result. As shown in the top right panel 
of Figure 8, the flux rope axis points predominantly in 
the GSE z-direction, with a small tilt in the ^-direction 
(as is clear from the direction of the k vector). This 
is consistent with the large tilt out-of-the-ecliptic found 
with the GCS fitting. For an ICME with such an orien- 
tation, it is expected that the cross-section in the ecliptic 
plane is relatively small and the cloud can be detected at 
Earth only if it propagates very close to the Sun-Earth 
line. The presence of a non-zero component of Bi inside 
the cloud indicates a non-zero impact parameter. Con- 
sidering the cloud is left-handed, this means that the 
center of the ejecta passed slightly west of the Sun-Earth 
line (meaning it passed on the — y side of Earth in GSE 
coordinates). The orientation of the shock normal (0.90, 
-0.40, -0.15) in GSE coordinates is also consistent with 
this picture. 



Mostl et al f 2009) was also performed, and it also yields 



a single, left-handed ejecta, with an axis making an an- 
gle of 10° with respect to the — z direction with a small 
impact parameter of 0.1 relative to the radial size of the 
flux rope. The cross-section map is shown in the bottom 
right panel of Figure 8. The view is approximately from 
ecliptic south with the Sun on the right-hand side of the 
picture. The Grad-Shafranov reconstruction technique 
has been shown to work best for magnetic clouds and 
magnetic cloud-like ejecta (as com pared to less regular 
ejecta, see Al-Haddad et al. (2012) for details) and this 
shows once again that the in situ measurements indicate 
a relatively typical flux rope-type ejecta, except for the 
higher plasma temperature. The axial field strength is 
found to be 13.5 nT, which is typical for a magnetic cloud 
or magnetic cloud-like ejecta at solar minimum. 

Returning to the left column of Figure 8, one may no- 
tice that before the ICME is crossed, the sheath density 
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Front 


Method 


Direction 


Arrival Time 


Speed at 1 AU 


HI2-A Track 1 


F$ 


5° 


05:00UT on 05/27 


630 km s _1 


HI2-A Track 1 


F<£ 


12° 


13:00UT on 05/27 


530 km s _1 


HI2-A Track 1 


HM 


24° 


09:00UT on 05/28 


370 km s -1 


HI2-A Track 2 


F<£ 


5° 


03:00UT on 05/28 


490 km s _1 


HI2-A Track 2 


HM 


5° 


11:00UT on 05/28 


390 km s _1 


HI2-B Track 


F<E> 


5° 


04:00UT on 05/28 


390 km s _1 


HI2-B Track 


HM 


24° 


09:00UT on 05/28 


360 km s _1 


HI1-A LEI 


F$ 




05:30UT on 05/28 


370 km s- 1 


HI1-A BC1 


F$ 




19:00UT on 05/28 


340 km s _1 


HI1-A LE2 


HM 


28° 


23:00UT on 05/28 


345 km s _1 


Shock 


Wind 




02:00UT on 05/28 


390 km s _1 


End of Sheath - Start Ejecta 


Wind 




20:00UT on 05/28 


390 km s _1 


Compression - End Ejecta 


Wind 




20:00UT on 05/29 


360 km s _1 



TABLE 1 

Predicted and measured arrival times and speed of the bright tracks observed by SECCHI on May 23-27, 2010 and in situ 
by Wind. LE refers to leading edge and BC to the beginning of the cavity. The predictions using HI1-A data assume a 

CONSTANT SPEED IN HI-2 (SEE SECTION 5.2 FOR DETAILS). 



decreases steadily while the magnetic field strength in- 
creases. These are features similar to those of a plasma 
depletion layer. A flow enhancement consistent with this 
expectation is shown in Figure 8, left column third panel, 
at the location of the second vertical guideline. In this 
interpretation, it is a draping effect. Draping around 
ICMEs o bserved at Eart h has been described previously 
(e.g., see |Liu et al.||2006| ). 

5.2. Prediction from HI Observations 

We use the measurements in HI2 combined with vari- 
ous estimates of the direction of the CMEs to calculate 
the arrival times at Earth of the bright fronts. The re- 
sults are summarized in Table 1. Using the F<3> approx- 
imation and a direction of 5° to 12° for the first track 
in HI2-A, the predicted arrival time is 12 to 20 hours 
earlier than any structure detected by Wind. This con- 
firms that the first track in HI2-A is not associated with 
CME1. In order to obtain a consistent arrival time and 
speed at 1 AU for HI2-A and HI2-B data, it is necessary 
to consider that the first track in HI2-A and the track 
observed in HI2-B correspond to CME2 with a direction 
no 24°. This substantiates the discussion in section 4.5, 
where we found that the first track observed in HI2-A 
corresponds to CME2. It should be noted that a CME 
passing 24° off the Sun-Earth line is not expected to re- 
sult in the type of clear flux rope measurements as seen 
in Figure 8, especially for a CME with such a large tilt 
with respect to the ecliptic. 

As explained in section 4.5 and further confirmed 
above, the second track in HI2-A corresponds to CME1. 
Assuming a direction of 5° for CME1 (due to the de- 
flection during the interaction, see section 6 for more 
details), the arrival time at Wind can be well reproduced 
using the F<I> approximation, but the predicted arrival 
speed is too high by more than 200 km s -1 . The op- 
posite is true with the HM approximation (arrival time 
too late by 7 hours but speed well reproduced). This 
might be because CME1 has a narrow, but non-negligible 
width, and, also because at large elongation angles, nei- 
ther method works well. 

We also use the observations in HI1 to predict the ar- 
rival times of CME1 and CME2 without taking into ac- 
count any additional acceleration or deceleration during 
the remaining of the CME-CME interaction. We simply 



assume that, past 0.35 AU, the CMEs propagate with 
constant speeds equal to the final speed determined from 
HI1 data. Results are also summarized in Table 1. As- 
suming a constant speed of 370 km s _1 , the leading edge 
of CME1 is predicted to arrive at 1 AU within four hours 
of the measured shock arrival time. The predicted arrival 
time for the magnetic ejecta is within two hours of the 
start of the magnetic cloud observed by Wind. 

Using the expansion rate observed before the CME- 
CME collision, CME1 is expected to arrive at 1 AU with 
a radius of about 0.11 AU, which is similar to the mea- 
sured radius of 0.1 AU. The same procedure yields an 
expected thickness of about 0.1 AU for the sheath in 
front of CME1, a slight underestimation of the measured 
sheath thickness. 

6. DISCUSSION: SCENARIO OF THE CME-CME 
INTERACTION 

Here, we attempt to explain simultaneously the in 
situ measurements showing a single ICME with a flux- 
rope type structure propagating close but to the west 
of the Sun-Earth line and the remote-sensing observa- 
tions showing the interaction of two CMEs. To do so, 
it is necessary to consider that, when the two CMEs in- 
teract, they are deflected away from each other, which 
ultimately results in the two CMEs separating from each 
other. 

As determined in sections 3.2 and 4.3, CME1 has an 
initial direction of about 12° and a half- width of about 
10°. CME2 has an initial direction of about 28° (see sec- 
tion 3.2). The collision between the two CMEs involves a 
relatively small cross-section (as sketched in the bottom 
panel of Figure 5). As the two CMEs collide, there is an 
eastward deflection of CME1 and westward deflection of 
CME2. CME1 is originally propagating along W12 and 
is observed at 1 AU with a very small impact parame- 
ter, consistent with a small western direction. Therefore, 
we can infer that the deflection of CME1 is of the order 
of 8-10°. The only numerical study of CME deflection 
during CME-CME interaction was performed by |Xiong| 
|et al.| ( [2009] ), focusing only on the latitudinal deflection. 
The authors found a maximum deflection when the two 
CMEs propagate 15° away from each other, as is the case 
here (see their Figure 12c). However, the deflection was 
found to be only 3-6° for their CME1. 

Initially, CME1 is compressed by the CME-CME inter- 
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action; because of the deflections, the two CMEs separate 
and CME1 is "free" to expand again. It has been pre- 
viously pro posed (see for example Figure 6 of Gulisano 
|et al.|2010| ) that a compressed magnetic cloud may over- 
expand to reach a size similar to what would be its size in 
the absence of the compression event. Compression may 
be due to a fast wind solar wind stream overtaking the 
magnetic cloud or an instance of CME-CME interaction. 
We believe that such an over-expansion is possible only 
if there is no other magnetic obstacle (another CME) at 
the back of CME1 to hinder it. Therefore, CME1 and 
CME2 must separate after the collision so that the mag- 
netic ejecta associated with CME1 over-expands to its 
expected size (that given by the relation found in sec- 
tion 4.3). The over-expansion results in CME1 having 
a typical radius at 1 AU and no in situ signatures of 
its compression, except, possibly, the higher proton tem- 
perature being present. This is also the reason why the 
measured width of CME1 at 1 AU is consistent with the 
prediction from HI1 data only. 

Due to its eastward deflection, Wind observes CME1, 
including its sheath and associated magnetic flux rope 
structure. CME2 is not observed in situ, having been de- 
flected further to the west. After the data gaps in SEC- 
CHI, CME1 is observed in HI2-A as the second track 
and in HI2-B as the first (and only) track. CME2 is 
observed in HI2-A as the first track, which shows a ge- 
ometrical acceleration consistent with a direction closer 
to STEREO- A than CME1. During the first 12 hours 
of propagation in HI2-A, there is only one bright feature 
observed in background-subtracted images. During this 
time, the signals from the two CMEs are superimposed 
and the track associated with CME2 "passes" that as- 
sociated with CME1. Because of the running-difference 
procedure, the tracks appear continuous in the classical 
J-map format, except for the "split" in the first track. 

7. CONCLUSIONS 

In this paper, we have analyzed in detail the remote- 
sensing observations of two successive CMEs in the in- 
ner heliosphere and the associated in situ measurements. 
The two CMEs were associated with eruptions from the 
same filament channel and they are clearly observed to 
interact in SECCHI/HIl-A images. The interaction is 
associated with a simultaneous decrease in speed of the 
second CME, and a compression in the radial direc- 
tion of the first CME. We have shown how, due to the 
two CMEs having different initial directions (by about 
10-15°), CME1 is deflected towards the Sun-Earth line 
(eastward) and CME2 away from Earth. At 1 AU, only 
CME1 is observed as an ICME with a radius of 0.11 AU 
and preceded by a sheath of thickness 0.15 AU. The 
ICME has the geometry of a flux rope with an axis along 
the —^-direction and all the characteristics of a magnetic 
cloud except an elevated plasma temperature. While the 
sheath has interesting features, including some evidence 
of draping around the ICME, there is no hint of a past 
CME-CME interaction in the in situ measurements, ex- 
cept perhaps for its high density and long (but not un- 
usual) duration. 

In our study of remote-sensing observations, we have 
made, or confirmed, a number of other findings. 1) We 
have determined that the radial expansion of the mag- 
netic ejecta with distance in the inner heliosphere can be 



approximated by a power law of exponent 0.82 (similar 
to th e studies by |Bothmer fc Schwenn||1994] |Liu et aL 



2005) and there is a weaker radial dependence for the 
sheath thickness. 2) We have found that the radius of 
the previously perturbed ICME at 1 AU is consistent 
with the expansion rate calculated in the inner helio- 
sphere between 0.1 and 0.35 AU (before the collision), 
confirming previous studies of compression of magnetic 
clouds by solar wind streams and other magnetic clouds 
( |Gulisano et al.||2010 ). 3) We have shown that each part 
of a CME observed in HI images can be assumed to prop- 
agate in the same direction, including CME tracks which 
are created through the running-difference procedure (for 
example, the "black" edge). 4) We have derived the tem- 
poral offset between the actual leading edge of a CME as 
seen in white-light images and the "black" edge as seen 
in running-difference images. This black edge is what is 
typically tracked in HI2 and used to estimate the CME 
arrival times from fitt ing t echniqu e s such as those from 



Rouillard et aL] ( |2008[ ) and |Lugaz| ( |2010 ). For the 2010 
May 23 CME, we found that the temporal offset increases 
with time and is approximately 2-3 hours in HI1 and 4-6 
hours in HI2. 

This study is also meant to demonstrate the difficul- 
ties of relying solely on remote-sensing observations to 
predict complex ICME events at Earth. While helio- 
spheric imagers have improved the predicti on of arrival 
times and arrival speed of CMEs (e.g., see |Davis et al.| 
2011||Liu et al.|2010b||Schreiner et al.|2012[ ), we are still 
struggling to understand complex events, such as the Au- 
gust 2010 series of CMEs or this series of two CMEs in 
May 2010. Here, for example, the tracks associated with 
the two CMEs do not appear to merge or cross with 
each other in the J-map made with STEREO-A data, 
whereas, in fact, the second, faster CME collided with 
the first, slower CME and may be ahead of it (but along 
a different direction) at 1 AU. In addition, there is only 
one bright track in HI2-B 

During the interaction, the speed of CME2 is found to 
decrease from close to 600 km s _1 to about 380 km s _1 . 
Comparing its speed in the HI1 FOV with the speed 
measured at 1 AU, it is clear that CME1 is only slightly 
accelerated during the collision. It would appear that 
the changes in the speeds of the two CMEs are consis- 
tent with a perfectly inelastic collision, but a more dedi- 
cated study is required to analyze in details the type of 
interaction. Such a study should take into consideration 
the expansion of the CMEs, their deflection, the internal 
magnetic field in the two CMEs, the reconnection be- 
tween the CMEs, and the possible presence of a shock 
ahead of CME2. Similarly to what was deter mined in 
the numerical simulation of |Lugaz et al.| ( |2005[ ), we find 
that, following the collision, the front of C1V1E2 propa- 
gates with a speed comparable with, but slightly higher 
than the speed of the front of CME1. This results in 
a contraction of CME1 after the collision, as its back 
moves faster than its front. We speculate that the speed 
of CME2 after the collision might be constrained by the 
speed of the leading edge of CME1 because the mag- 
netic tension and pressure inside CME1 limit the rate of 
contraction of CME1. In this interpretation, CME1 con- 
traction speed constrains the speed of the front of CME2 
to be close to, but slightly higher than the speed of the 
front of CME1. 
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In the future, it will be important to perform more 
dedicated numerical studies of the interaction of two (or 
more) CMEs with different orientations, sizes and longi- 
tudinal (and latitudinal) separations in three dimensions, 
with high resolution and including syntheti c HI images, 
a com bination of the detailed si mulation of |Lugaz et aL 
(120051) with the cases studied inlXiong et al.l (12009). 



The authors would like to thank the reviewer for 
his/her useful comments which have helped make our 
manuscript clearer. The research for this manuscript 
was supported by the following grants: NASA NAS5- 
00132, NNX10AQ29G, and NNX12AB28G and NSF 



AGS-1239699, and AGS-1239704. I. R. would like to 
acknowledge the support from CAS grant 2011T2J01 at 
the YNAO and NSF grant AGS-0639335 (CAREER) at 
the IfA. M. T. acknowledges the Austrian Science Fund 
(FWF): V195-N16. C. M. acknowledges support from 
a Marie Curie International Outgoing Fellowship within 
EU FP7 and from the EU FP7 program under grant 
agreement n° 263252 [COMESEP]. SoHO and STEREO 
are projects of international cooperation between ESA 
and NASA. The SECCHI data are produced by an inter- 
national consortium of NRL, LMSAL, and NASA GSFC 
(USA), RAL, and U. Birmingham (UK), MPS (Ger- 
many), CSL (Belgium), IOTA, and IAS (France). 



REFERENCES 



Abraham-Shrauner, B. 1972, J. Geophys. Res., 77, 736 
Al-Haddad, N., Nieves-Chinchilla, T., Mostl, C, et al. 2012, Solar 

Phys., in revision 
Bothmer, V., & Schwenn, R. 1994, Space Sci. Rev., 70, 215 
Burlaga, L., Berdichevsky, D., Gopalswamy, N., Lepping, R., Sz 

Zurbuchen, T. 2003, J. Geophys. Res., 108, 2 
Burlaga, L., Sittler, E., Mariani, F., & Schwenn, R. 1981, J. 

Geophys. Res., 86, 6673 
Burlaga, L. F., Behannon, K. W., & Klein, L. W. 1987, J. 

Geophys. Res., 92, 5725 
Burlaga, L. F., Plunkett, S. P., & St. Cyr, O. C. 2002, J. 

Geophys. Res., 107, 1 
Colaninno, R. C., & Vourlidas, A. 2009, Astrophys. J., 698, 852 
Davies, J. A., Harrison, R. A., Rouillard, A. P., et al. 2009, 

Geophys. Res. Lett., 36, L02102 
Davies, J. A., Harrison, R. A., Perry, C. H., et al. 2012, 

Astrophys. J., 750, 23 
Davis, C. J., de Koning, C. A., Davies, J. A., et al. 2011, Space 

Weather, 90, S01005 
de Koning, C. A., Pizzo, V. J., & Biesecker, D. A. 2009, Solar 

Phys., 256, 167 
Eyles, C. J., Harrison, R. A., Davis, C. J., et al. 2009, Solar 

Phys., 254, 387 

Farrugia, C. J., Jordanova, V. K., TThomsen, M. F., et al. 2006a, 

J. Geophys. Res., Ill 
Farrugia, C. J., Matsui, H., Kucharek, H., et al. 2006b, Advances 

in Space Research, 38, 498 
Goldstein, H. 1983, in NASA Conference Publication, Vol. 228, 

NASA Conference Publication, 731-733 
Gopalswamy, N., Yashiro, S., Kaiser, M. L., Howard, R. A., &; 

Bougeret, J.-L. 2001, Astrophys. Journ. Lett., 548, L91 
Gulisano, A. M., Demoulin, P., Dasso, S., Ruiz, M. E., & Marsch, 

E. 2010, Astron. Astrophys., 509, A39 
Harrison, R. A., Davis, C. J., Eyles, C. J., et al. 2008, Solar 

Phys., 247, 171 

Harrison, R. A., Davies, J. A., Rouillard, A. P., et al. 2009, Solar 
Phys., 256, 219 

Harrison, R. A., Davies, J. A., Mostl, C, et al. 2012, Astrophys. 
J., 750, 45 

Howard, R. A., Moses, J. D., Vourlidas, A., et al. 2008, Space Sci. 
Rev., 136, 67 

Howard, T. A., & DeForest, C. E. 2012, Astrophys. J., 746, 64 
Howard, T. A., & Tappin, S. J. 2009, Space Sci. Rev., 147, 31 
Hu, Q., & Sonnerup, B. U. 6. 2001, Geophys. Res. Lett., 28, 467 
Hundhausen, A. J. 1993, J. Geophys. Res., 98, 13177 
Kahler, S. W., & Webb, D. F. 2007, J. Geophys. Res., 112, 9103 
Kaiser, M. L., Kucera, T. A., Davila, J. M., et al. 2008, Space Sci. 
Rev., 136, 5 

Kilpua, E. K. J., Pomoell, J., Vourlidas, A., et al. 2009, Annales 

Geophysicae, 27, 4491 
Lepping, R. P., Acuna, M. H., Burlaga, L. F., & M., F. W. 1995, 

Space Sci. Rev., 71, 207 
Lin, R. P., Anderson, K. A., Ashford, S., et al. 1995, Space Sci. 

Rev., 71, 125 

Liu, Y., Davies, J. A., Luhmann, J. G., et al. 2010a, Astrophys. 

Journ. Lett., 710, L82 
Liu, Y., Richardson, J. D., & Belcher, J. W. 2005, Planet. Space 

Sci., 53, 3 



Liu, Y., Richardson, J. D., Belcher, J. W., Kasper, J. C, &; 

Skoug, R. M. 2006, J. Geophys. Res., Ill, 9108 
Liu, Y., Thernisien, A., Luhmann, J. G., et al. 2010b, Astrophys. 

J., 722, 1762 

Liu, Y. D., Luhmann, J. G., Mostl, C, et al. 2012, Astrophys. 

Journ. Lett., 746, L15 
Lopez, R. E. 1987, J. Geophys. Res., 92, 11189 
Lugaz, N. 2010, Solar Phys., 267, 411 

Lugaz, N., Hernandez-Charpak, J. N., Roussev, I. I., et al. 2010, 

Astrophys. J., 715, 493 
Lugaz, N., Kintner, P., Mostl, C, et al. 2012, Sol. Phys., 279, 497 
Lugaz, N., Manchester, W. B., & Gombosi, T. I. 2005, Astrophys. 

J., 634, 651 

Lugaz, N., Manchester, W. B., Roussev, I. I., & Gombosi, T. I. 

2008a, J. Atmos. Solar-Terr. Phys., 70, 598 
Lugaz, N., Manchester, W. B., Roussev, I. I., Toth, G., & 

Gombosi, T. I. 2007, Astrophys. J., 659, 788 
Lugaz, N., Vourlidas, A., & Roussev, I. I. 2009a, Annales 

Geophysicae, 27, 3479 
Lugaz, N., Vourlidas, A., Roussev, I. I., et al. 2008b, Astrophys. 

Journ. Lett., 684, Llll 
Lugaz, N., Vourlidas, A., Roussev, I. I., & Morgan, H. 2009b, 

Solar Phys., 256, 269 
Manchester, W. B., Vourlidas, A., Toth, G., et al. 2008, 

Astrophys. J., 684, 1448 
Mostl, C, Farrugia, C. J., Biernat, H. K., et al. 2009, Solar Phys., 

256, 427 

Mostl, C, Temmer, M., Rollett, T., et al. 2010, Geophys. Res. 

Lett., 37, L24103 
Mostl, C, Rollett, T., Lugaz, N., et al. 2011, Astrophys. J., 741, 

34 

Odstrcil, D., & Pizzo, V. J. 2009, Solar Phys., 259, 297 
Odstrcil, D., Vandas, M., Pizzo, V. J., & MacNeice, P. 2003, in 
American Institute of Physics Conference Series, Vol. 679, Solar 
Wind Ten, ed. M. Velli, R. Bruno, F. Malara, & B. Bucci, 
699-702 

Robbrecht, E., & Berghmans, D. 2004, Astron. Astrophys., 425, 
1097 

Rollett, T., Mostl, C, Temmer, M., et al. 2012, Solar Phys., 276, 
293 

Rouillard, A. P. 2010, J. Atmos. Solar-Terr. Phys. 
Rouillard, A. P., Davies, J. A., Forsyth, R. J., et al. 2008, 

Geophys. Res. Lett., 35, L10110 
Savani, N. P., Davies, J. A., Davis, C. J., et al. . 2012, Solar 

Phys., 279, 517 

Schmidt, J., & Cargill, P. 2004, Annales Geophysicae, 22, 2245 
Schreiner, S., Cattell, C, Kertsen, K., & Hupach, A. 2012, Solar 

Phys., doi:10.1007/sll207-012-9936-5 
Sheeley, N. R., Walters, J. H., Wang, Y.-M., & Howard, R. A. 

1999, J. Geophys. Res., 104, 24739 
Shen, C, Wang, Y., Wang, S., et al. 2012, Nature Phys., in press 
Song, H.-Q., Chen, Y., Li, G., Kong, X.-L., and Feng, S.-W.: 

2012, Physical Review X 2, 021015. 
Sonnerup, B. U. O., & Cahill, Jr., L. J. 1967, J. Geophys. Res., 

72, 171 

Temmer, M., Vrsnak, B., Rollett, T., et al. 2012, Astrophys. J., 
749, 57 



14 



LUGAZ ET AL. 



Thernisien, A., Vourlidas, A., & Howard, R. A. 2009, Solar Phys., 
256, 111 

Thernisien, A. 2011, Astrophys. J. Supp., 194, 33 
Vandas, M., Fischer, S., Dryer, M., et al. 1997, J. Geophys. Res., 
102, 22295 

Wang, Y. M., Wang, S., & Ye, P. Z. 2002, Solar Phys., 211, 333 
Wang, Y. M., Ye, P. Z., Wang, S., & Xue, X. H. 2003, Geophys. 

Res. Lett., 30, 33 
Webb, D. F., Howard, T. A., Fry, C. D., et al. 2009, Solar Phys., 

256, 239 

Wood, B. E., Howard, R. A., Phmkett, S. P., &; Socker, D. G. 

2009, Astrophys. J., 694, 707 
Wu, C.-C., Fry, C. D., Wu, S. T., Dryer, M., & Liou, K. 2007, J. 

Geophys. Res., 112, A09104 



Wu, S. T., Wang, A. H., & Gopalswamy, N. 2002, in ESA SP-505: 

SOLMAG 2002. Proceedings of the Magnetic Coupling of the 

Solar Atmosphere Euroconference, 227-230 
Xie, H., Gopalswamy, N., Manoharan, P. K., et al. 2006, J. 

Geophys. Res., Ill, 1103 
Xiong, M., Zheng, H., & Wang, S. 2009, J. Geophys. Res., 114, 

A11101 

Xiong, M., Zheng, H., Wang, Y., & Wang, S. 2006, J. Geophys. 
Res., Ill, A08105 



